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Abstract 

We study the behavior of a point particle incident from the left on a slab of 
a randomly diluted triangular array of circular scatterers. Various scattering 
properties, such as the reflection and transmission probabilities and the scat- 
tering time are studied as a function of thickness and dilution. We show that 
a diffusion model satisfactorily describes the mentioned scattering properties. 
We also show how some of these quantities can be evaluated exactly and their 
agreement with numerical experiments. Our results exhibit the dependence 
of these scattering data on the mean free path. This dependence again shows 
excellent agreement with the predictions of a Brownian motion model. 

PACS number(s): 05.40. +j, 05.45.+b, 05.60. +w 

I. INTRODUCTION 

The Lorentz gas |3J is one of the fundamental kinetic theory models || , and many of its 
ergodic properties are well known . In a previous paper [Q we investigated the scattering 
and transmission properties of an array of disks centered on a finite triangular lattice, a 
Lorentz slab. There, we studied the validity of approximating the motion of such a particle 
by a diffusive process ||. Specifically, we considered samples of finite thickness L along 
the x axis, but infinite along the y axis. We then looked at the behavior of a particle 
incident on the sample from the left and measured the probability T of being transmitted 
to the right, the reflection probability R and the average residence time (r). We found a 
very satisfactory agreement with the predictions of one dimensional Brownian motion for all 
these quantities when the distance between scatterers is small enough to prevent any particle 
from moving in an arbitrarily long straight line (the so-called finite horizon case). When 
such unbounded motion is possible we found that normal diffusive behavior did not arise. 
Rather, a complex pattern of logarithmic corrections was found for the various quantities 
of interest. A modification of the diffusive model considering Levy walks as suggested in 
Refs. 0,0] explains such features. 
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Here we study the situation in which the system is randomly diluted, but the scatterers 
are still placed on the the sites of a triangular lattice. That is, we consider the case in which a 
fraction / of the cells of the periodic array is occupied by scatterers and the rest is empty. In 
such a system, the particle can always take steps of arbitrary length, but the horizon is said 
to be finite if the distribution of free paths has an exponential cutoff. This happens when 
the corresponding system with / = 1 has finite horizon, i.e., the only large steps that occur 
in the diluted system are related to the absence of a large number of scatterers, which is an 
exponentially improbable event. In this paper, we restrict ourselves to this finite horizon 
case. 

The importance of looking at the randomly diluted case, is that it provides a controlled 
way of varying the mean free path — and hence the diffusion constant — over an arbitrarily 
large range, without affecting the property of having a finite horizon and hence a normal 
diffusive behavior. This therefore allows us to check the correspondence with diffusion 
models much more thoroughly. Further, it shows that the correspondence between classical 
deterministic motion and Brownian motion, as proved rigorously in Ref. || for the case of 
periodic billiards, does not, in fact, require strict periodicity. 

This paper is organized as follows: In Sec. |I| we describe in detail how the numerical 
experiments are performed. In particular, we define two different ways in which we introduce 
disorder in the system, namely quenched and annealed. In Section [TTT| , we discuss the average 
transport properties, such as transmission probability, mean free path and average scattering 
time. We find that the last two quantities can be evaluated exactly using a relationship due 
to Kac ||. We display the numerical results for these quantities as well as the predictions 
using the Kac formula and the diffusion model. For the sake of completeness, we reproduce 
the derivation of the Kac formula in Appendix A. In Sec. |V| we discuss the distribution 
functions of free paths, of residence times and of heights of exit of transmitted particles. In 
Appendix B we derive the latter two distributions for the diffusion model. Finally in Sec. [V| 
we present our conclusions. 

II. MODEL SETUP 

The geometric arrangement of the scatterers in the undiluted (/ = 1) system is the 
following: each scatterer is a disk of unit radius, the centers of which form a triangular 
lattice. The slab is infinite in the vertical direction and is characterized by the number L 
of columns and the minimal separation w between the disks. Dilution is then obtained by 
occupying only a fraction / of the sites of the array with disks of unit radius. A typical 
scatterer configuration is shown in Fig. |I[ The left and right "sawtooth" borders are the 
outer sides of the hexagonal cells attached to each site of the lattice. 

Particles are launched from randomly chosen positions along the left border. Each inci- 
dent particle has a different impact parameter b, defined here as the distance between the 
initial position and the horizontal line passing through the center of the scatterer in the 
cell. The angles of incidence 9 measured with respect to the side from which the particle is 
launched, are distributed in the interval [0, 7r] in such a way as to make cos 9 uniformly dis- 
tributed. This choice reproduces initially the Liouville measure in the Birkhoff coordinates. 
The particles move freely except for elastic collisions at the boundary of the disks. 
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In the undiluted case (/ = 1), if the separation between scatterers w is small, < w < 
w c = (4/ a/3 — 2) = 0.3094 . . ., the length of free motion of the particles is bounded and they 
"see" a finite horizon. All the numerical experiments discussed in this paper are performed 
in this range. On the other hand, in the diluted case, arbitrarily long paths without collisions 
can exist also for < w < w c , but their contribution to the diffusion constant remains finite. 

As far as the dilution process itself is concerned, we realize it in two different ways cor- 
responding to the usual distinction between quenched and annealed disorder. The annealed 
disorder is obtained by choosing with probability / the cell to be occupied at the moment 
at which the particle enters the cell. Thus, when the particle eventually revisits a given cell, 
its occupancy status may be different. This way of introducing disorder involves a simulta- 
neous average over the dynamical and the disorder variables. On the other hand, we also 
performed dilution in the more realistic quenched case, in which a sample is first created, 
for which all cells are either occupied or empty with probability / and averages are taken 
over many realizations of disorder. 



III. AVERAGE TRANSPORT PROPERTIES 

We begin this Section by showing a derivation of a formula for the mean free path A as 
a function of the geometrical parameter w and dilution /. This quantity can be evaluated 
in terms of phase space integrals using Kac formula ||, which assumes ergodicity. Indeed, 
as shown in Appendix A, for a particle moving freely at unit velocity in a bounded domain 
of area S, the average return time (T) to a boundary segment of length P is 

(T) = (1) 

For the Lorentz gas, choosing the boundary segment as the perimeter of all disks, the average 
return time coincides with the mean free time between two collisions. Since the particles 
move with unit velocity, this also coincides with the mean free path A. Considering a finite 
array of Wigner-Seitz hexagonal cells with L columns and M rows, the total area of the 
domain is LMC, where C = (2 + w) 2 \f?>/2 is the cell area. From this we must subtract the 
area occupied by the unit radius disks, i.e. nfLM. The total perimeter is P = 2nfLM. 
Hence, from Eq. ([!]) and the above considerations 



A = ^(- + 2) 2 -^. (2) 



Note that this holds for any value of M and hence extends trivially to the case we consider, 
where M is infinite. The usual derivation of Eq. (|l|) only applies to the quenched case, 
i.e. to a scatterer configuration that is fixed in time. However, the arguments we give in 
Appendix A show how to extend it to the annealed well. Numerical experiments 

confirm this finding, for both quenched and annealed disorder. In Fig. [2] we display the 
results for quenched disorder, since the annealed data are identical. 

The slab is infinite in the y direction and the collisions are elastic, hence every particle 
that enters the slab must eventually leave it, except for a set of zero measure which goes 
asymptotically to bounded orbits inside the slab. Thus, in practice, a particle that enters 
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the slab collides with some of the obstacles and is ultimately transmitted or reflected. From 
Eq. (jl]) one can obtain the average residence time (r) as a function of w, f and L, if one now 
takes the border of the slab as the boundary segment. For a slab with L columns and M 
rows the perimeter is M4(2 + w)f \/3 and the total area of the domain is as before. Hence 

This equation defines the quantity B. As above (see Eq. (|2|)), the result is independent of 
M. We verified numerically the linear growth of (r) with L and found it to be valid, as 
expected, also for small values of L. From these experiments we obtained the values of B 
shown in Fig. ^ for quenched disorder together with Eq. . The agreement in the annealed 
case (not shown) is even more satisfactory. Note that B does not diverge as / — > 0, as 
opposed to the mean free path. 

Let us now compare these exact results with the predictions of a model based on Brownian 
motion. We assume that the particles incident from the left penetrate a given distance a 
into the sample, after which they forget everything about the way they were injected into 
the system and diffuse with diffusion constant D. Under these assumptions, the quantities 
of interest can be evaluated exactly in terms of a and D and compared with the exact results 
obtained above. Thus, the average time to reach either side starting at a distance a from 



the left side is given by [110,11 



(r)(I) = ^, (4) 

where I = v3(l + w/2) is the horizontal separation between columns of the slab. One sees 
therefore that the form of the L dependence in the diffusive case is slightly different from the 
exact one, Eq. (HD, since in that case the proportionality to L is exact over the whole range 
of L. On the other hand, in Eq. (|4]), a constant term appears which is negligible in the limit 
aft < L. Since this is the limit for which diffusion is expected to be a valid description, this 
is not a serious problem. On the other hand, this allows to evaluate the ratio a/D exactly 
in terms of / and w, via 

a 2B(f,w) 

D = — — (5) 

It is well known that it is not possible to obtain exact evaluations of the same sort for the 
diffusion constant D, since its value can be shown to depend in a detailed manner on the 
specific dynamics involved. Indeed, 



D = lim J x(0)x(t)e- et dt, (6) 

where the bar denotes an average over the realizations of the disorder. From this it follows 
that an analytical expression for D in terms of simple phase space expressions is impossible. 
However, an order of magnitude estimate for D can be given. Since the mean square distance 
grows linearly with time, and since the only microscopic length scale is the mean free path, 
we are led to estimate Dt by the square of the mean free path, where t is the time needed 
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for the particle to cover a mean free path. In our system, the velocity of the particle is 
constant and equal to one, therefore the result is that D is of the same order as a, both 
being of the order of the mean free path. This is indeed consistent with Eq. (Bl). It should 
be noted that the ratio of the two dynamics dependent quantities a and D depends only on 
the geometric features of the system. 

Let us now turn to another average transport property, for which no exact expression 
is available, namely the transmission probability T as a function of L. One finds that 
asymptotically, in the limit L ^> 1, 

T(f,w,L) = — - — (7) 

In Fig. [|, we show that A grows almost linearly with A, for large values of A. This can 
readily be understood in terms of the diffusive model. Indeed, in the case of a diffusing 
particle being injected at a distance a from the left-hand side of the slab, the probability 
that it be transmitted to the right-hand side without being absorbed first at the left-hand 
side T d is well-known to be iTUIjn 



UL) = ±. (8) 

From this we see that A can be identified with a quantity which, as argued above, is expected 
to scale as the mean free path. Note further that whereas (r) and A are identical in the 
quenched and annealed cases, this is not true for the transmission probability, which shows 
significant differences for the two forms of disorder. 



IV. DISTRIBUTION FUNCTIONS 

We have also studied the distributions of residence times of transmitted (respectively 
reflected) particles, of free paths and of the heights at which transmitted (respectively re- 
flected) particles leave the system. 

Of these, only the distribution of free paths is accessible to exact theoretical treatment. 
Indeed, this is nothing but an equilibrium property, and can in principle be evaluated using 
an integral over the perimeter of the billiard with the Birkhoff measure. However, this is 
not really practicable, so we have not attempted it. The distribution of mean free paths 
is quite irregular (see Fig. |5|). This irregularity arises from the presence of resonances and 
the fact that there are forbidden distances. Note that, as remarked in the Introduction, the 
distribution decays exponentially since the horizon is finite. From the above characterization 
of the mean free path distribution as an equilibrium property, it readily follows that they 
are independent of whether the average is taken as quenched or annealed. 

Now let us consider the distribution of residence times, which is not susceptible of an 
exact evaluation. However, it can be computed in the diffusion model. It is clearly sufficient 
to compute the distribution of residence times for transmitted particles ptrans upon starting 
at a, since the corresponding distribution for reflected particles p re fi arises by substituting 
a by £L — a in p tra ns- Defining the scaled variables r = Dr / (£L) 2 and a = a/(£L), the final 



result is, as shown in Appendix B [11 
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Ptrans{r; a) = 2tt n sin 7m(l - a) e 71 nr , (g) 



n=l 



Integrating, one obtains Eq. @ for the average time a diffusing particle takes to reach either 
side. We fit this to the numerical data as follows: one first considers the average time that 
a transmitted particle remains in the sample. This is obtainable from the experimental 
distribution on the one hand, but is also given by 

i 2 L 2 

(r)trans = ■ (10) 

This allows therefore a determination of D from the data, which for / = 1 is in good 
agreement with published data ||. From this and Eq. @ together with Eq. @, one also 
obtains a value of the parameter a. One is then in a position to plot the theoretical curve 
Eq. (|^) together with the empirical data. This is shown in Fig. ^|, for an occupation fraction 
/ = 0.5. The agreement is quite good. The apparent shift between the theoretical curve 
and the data in reflection can be traced back to an issue of normalization involving particles 
reflected after a very few bounces, which therefore do not show diffusive behavior. Thus, 
the empirical distribution has a short-time cut-off at larger times than the corresponding 
diffusive model. 

Finally, we measure the distribution of heights of transmitted particles. This distribution 
can be computed in the diffusive model. The result for transmitted particles is, as shown in 
Appendix B, 

p( v ) = Sin7r(1 ~ a) I , (11) 

2a cosily — cos7r(l — a) 

where rj is the scaled height yj [IV) and a is as above. In Fig. [7| we show p(rj) for transmitted 
particles. The agreement is excellent. We do not show the behavior for reflected particles 
as it is dominated by rapid reflections. 



V. CONCLUSIONS 

Summarizing, we have studied transport properties of finite size samples of Lorentz gases 
in a situation in which the mean free path can be varied over a large range of values without 
affecting the property of normal diffusion. The mean free path was varied over a considerable 
range (slightly more than an order of magnitude), and the agreement with the Brownian 
motion model was satisfactory throughout. We have also shown how the diffusion constant 
and the penetration depth, which were the two dynamical parameters of our model, vary 
with dilution, and hence with the mean free path. 
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APPENDIX A: 

In this Appendix, we show in detail how the exact relations Eq. (Q) and Eq. @ are 
obtained. To this end, we first recall the derivation of a general formula due to Kac. Consider 
a (2N — l)-dimensional energy shell in a 2A^-dimensional Hamiltonian phase-space and select 
a (2N — 2)-dimensional Poincare surface that intersects all (or nearly all) trajectories on the 
energy shell. If such a surface cannot be found, then the phase space integrals below must 
be restricted to that part of the phase space which can be reached from the surface. 

Parametrize each point (p, q) on the energy surface using the last point (p, q) on the 
Poincare surface that lies on the trajectory passing through (p,q). Denote by Te(p, q) the 
time necessary to reach (p,q) starting from (p,q). This defines the canonical coordinate 
transformation 

d N pd N q = d N - 1 pd N ~ 1 qdTdE. (Al) 
The constant energy Eq volume is hence 

J d N pd N q6[E - H(p,q)} = J d N ~ 1 p d N ' x q T Eo (p, q) . (A2) 

Denoting by A/", the total (2N — 2)-dimensional phase space volume of the Poincare surface, 
one immediately obtains the Kac formula for the average time to return to the surface: 

(Teo) = AT 1 J 5[E - H(p, q)}d N pd N q, (A3) 

N= f d N ~ l pd N - l q. 
Je=e 

To evaluate these integrals in the case of billiards for = 2 (of which the Lorentz gas with 
quenched disorder is a particular case), we take the Hamiltonian to be p 2 /2 and Eq to be 
1/2. From this follows 

J d 2 p d 2 q5[l/2 - p 2 /2] = 2tt S, (A4) 

where S is the area of the billiard. As a Poincare surface we introduce an arbitrary subset of 
the billiard boundary having perimeter P, with the usual Birkhoff coordinates as variables 
p and q. One then finds 

/ dpdq = P cos 9 d9 = 2P, (A5) 

JE=E J-tv/2 

From these two equations one derives Eq. ([TJ) 
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The above remarks are clearly limited to the case of quenched disorder. The general- 
ization to the annealed case can be made as follows: to every point (p,q) on the Poincare 
surface, add a doubly infinite sequence of zeroes and ones (cfe)fcL_ oo) which we denote by a, 
with all (jfc independently distributed and equal to one with probability /. We now define the 
dynamics as follows: The orbit starts from (p, q) and the cell at t = is occupied or empty 
according to the value of do. The orbit then proceeds until it leaves the cell. The status of 
the next cell is then decided according to the value of o"i and so on. Although discontinuities 
arise when a trajectory crosses a vertex of the fundamental cell, the dynamics is still given 
by a canonical map. Thus, for a fixed sequence a the formula Eq. ( |A3| ) applies. However, 
for a fixed sequence, we cannot easily compute the phase space volumes involved. Since 
we are only interested in the average of (T) (Zf) over all values of a, we note the following: 
the denominator in Eq. ( |A3| ) is independent of a, so that it is sufficient to average over the 
numerator. This average can be performed separately for each cell, and the result therefore 
follows trivially. 



APPENDIX B: 

In this Appendix, we derive Eq. ([|) and Eq. (|TT|). To this end, we first derive an analogue 
of Eq. (|9D for the probability that a particle first exits the slab on the right side at a height 
y and at time t. In order to simplify the notation, we first go over to scaled variables 
£ = x/(£L), r] = y/(£L), a = a/(£L) as well as r = Dt/(£L) 2 . All distribution functions are 
further rescaled in such a way as to remain normalized. In these variables, this probability 
is given by the following expression: 



PtransiV, T ) = ~ 



(Bl) 



where Po(£,,v) is the solution of the following problem: 

P (tr,;0)=S((-a)S( V ) P„(0,t,;t) = P (1,v,t) = 0. (B2) 

This is now solved by developing the delta function in eigenmodes of the Laplacian satisfying 
the boundary conditions in Eq. (|B2|), that is 



00 r°° rlk r i 

P (^ri;T) = 2 Y sinn7r£sinmra / — e ife? exp -(ttV + k 2 )r] . (B3) 

n=l 7-oo27T 



The resulting expression for ptransiv t) can now be integrated either over 77 or over r to 
yield Eq. (|j|) or Eq. (11) respectively. In either case, obtaining the results in the text is now 



a matter of straightforward algebra. 
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FIGURES 



FIG. 1. Slab of diluted scatterers in a triangular array with L = 10, w = 0.2, and / = 0.7. 
Point particles enter the slab on the left side saw tooth border. 



FIG. 2. Mean free path shifted by tt/2 as a function of the dilution / for quenched disorder. 
The curves are the r.h.s. of Eq. (^) again shifted by ir/2. Each point represents the average over 
10 8 collisions. 



FIG. 3. Dependence of B on / for quenched disorder. The curves are Eq (||). The values of B 
were found by adjusting Eq. (||) to the experimental data obtained by letting 10 7 particles travel 
through slabs of lengths going from 100 to 1500. 



FIG. 4. Dependence of A on the mean free path A for w = 0.2. Both quenched and annealed 
disorder are shown. The values of A were found by adjusting Eq. (|^) to the experimental data 
obtained by letting 10 particles travel through slabs of lengths going from 100 to 1500. 



FIG. 5. Free path distribution for w = 0.2, / = 0.5 for both quenched and annealed disorder. 
The curves are practically superimposed. The distributions were found from 10 8 collisions in each 
case. 



FIG. 6. Distribution of residence times for reflected (a) and transmitted (b) particles for 
quenched disorder, w = 0.2 and / = 0.5. The continuous curve represents the fit described in the 
text. The experimental distributions were found by letting 10 7 particles travel through a slab with 
100 columns. 



FIG. 7. Height distribution of transmitted particles for quenched disorder, w = 0.2 and 
/ = 0.2. The continuous curve is the plot of Eq. (11). The experimental distribution was found by 
letting 10 7 particles travel through a slab with 100 columns. 
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